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We report Suzaku observations of the galaxy cluster Abell 1795 that extend to r2oo ~ 2 Mpc, the radius 
within which the mean cluster mass density is 200 times the cosmic critical density. These observations 
are the first to probe the state of the intraclustcr medium in this object at r > 1.3 Mpc. We sample 
two disjoint sectors in the cluster outskirts (1.3 < r < 1.9 Mpc) and detect X-ray emission in only one of 
. them to a limiting (3 a) soft X-ray surface brightness of -Bo. 5-2 koV = 1-8 X 10~ 12 erg s _1 cm" 2 deg -2 , a 

C*~) ' level less than 20% of the cosmic X-ray background brightness. We trace the run of temperature with 

\q ■ radius at r > 0.4 Mpc and find that it falls relatively rapidly (Td op rojectcd °c r~°- 9 ), reaching a value about 

C^) ' one third of its peak at the largest radius we can measure it. Assuming the intracluster medium is in 

0^ I hydrostatic equilibrium and is polytropic, we find a polytropic index of T = 1.3l ; 2 and we estimate a mass 

of 4.11q 3 x 10 14 M Q within 1.3 Mpc, somewhat (2.7 a) lower than that reported by previous observers. 
However, our observations provide evidence for departure from hydrostatic equilibrium at radii as small as 
r ~ 1.3 Mpc « rsoo in this apparently regular and symmetrical cluster. 
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1. Introduction the vicinity of r2oo is difficult to observe efficiently with 

Chandra and XMM-Newton because of their relatively 
Spatially resolved X-ray spectroscopy is essential for high (and, in the latter case, time-variable) levels of in- 
mapping the distribution of mass, intracluster plasma, strumental background. 

and heavy elements in galaxy clusters. Although Chandra The study of cluster outskirts is important for several 

and XMM-Newton have provided a wealth of information reasons. First, a direct X-ray measurement of total clus- 

011 these quantities in the central regions of clusters, rel- ter mass (e.g., within r2oo) requires knowledge of the ICM 

atively little is known about their outer regions, at radii temperature profile at this radius. This is especially im- 

greater than ~ 0.5r2oo- Here r2oo is the radius within portant to facilitate comparison between X-ray and weak- 

which the mean cluster density is 200 times the critical lensing mass profiles, as the latter are difficult to measure 

cosmic density. In the simplest case of the collapse of in the inner regions of clusters. Second, structure and 

a spherical density perturbation in an Einstein-deSittcr cluster formation models make definite predictions about 

universe, it is expected that the dark and bayronic mat- total mass and ICM distributions to these radii (Navarro, 

ter within r2oo will have reached dynamical equilibrium Frenk & White 1997; Borgani et al. 2004; Roncarclli ct al. 

(Evrard, Metzler & Navarro 1996). Cluster emission from 2006), so measurements in this radial range offer a test of 
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these models. Finally, to the extent that V2oo marks the 
approximate boundary of the dynamically relaxed region 
of a cluster, one might expect to find inhomogencitics in 
the ICM that could provide a direct view of the accretion 
processes by which clusters grow. 

ROSAT has been used to characterize the X-ray surface 
brightness of the outer regions (r ~ r2oo) of individual clus- 
ters (Vikhlinin, Forman & Jones 1999) and, using stack- 
ing techniques, of ensembles of objects (Neumann 2005). 
Chandra surface brightness profiles in the range rsoo < 
t < ?'20o for a sample of relatively distant clusters have 
been studied by Ettori & Balestra (2008). Recent Suzaku 
observations have successfully probed out to and beyond 
r2oo for a handful of clusters (George et al. 2008; Fujita 
et al. 2008; Reiprich et al. 2008). 

In this paper we present Suzaku measurements of 
the surface brightness and temperature of the intraclus- 
ter medium of the X-ray luminous (Lx.boi = 1.4 x 10 45 
erg s _1 ) cluster Abell 1795. This is a relatively hot 
(spatially-averaged kT = 5.3 keV) cluster with a cool core 
(Horner 2001; Vikhlinin et al. 2006; Snowden et al. 2008). 
Although Chandra shows a 40-arcsec-long (50 /i^, 1 kpc) 
plume "trailing" the brightest cluster galaxy (Fabian et al. 
2001), outside of the very center the cluster is quite regu- 
lar and appears relaxed. From the observed temperature, 
and scaling relations derived from simulations (Evrard, 
Metzler & Navarro 1996; Arnaud & Evrard 1999) we es- 
timate rsoo = 1.3 h^Q Mpc, consistent with the Chandra 
measurements of Vikhlinin et al. (2006), and r2oo = 
Mpc (26') for our assumed cosmology 1 . We exploit the 
relatively low and stable Suzaku background to trace the 
temperature and ICM density to the vicinity of f2oo- 

As noted by Bryan & Norman (1998), if the virial radius 
fviriai = J"a„ is derived from the collapse of a spherical 
top-hat density perturbation, assuming that the cluster 
has just virialized, then for an Einstein-deSitter cosmology 
A„ = 18-7T 2 w 180, independent of redshift. In general, 
however, A„ is function of cosmology and redshift, and 
for our ACDM cosmology at the redshift of Abell 1795 
A v m 100, so the virial radius is considerably larger than 
r 2 oo- For Abell 1795 we expect ri o ~ 1.35r 2 oo- Here A„ 
is the ratio of the mean cluster mass density within ta„ 
to the critical cosmic density. 

2. Observations, Data Reduction &; Analysis 

2.1. Observations 

Abell 1795 was observed on 2005 Dec 10, during the 
Suzaku performance verification phase, with a series of 
5 overlapping pointings, illustrated in Figure 1. The re- 
sults reported here were obtained with XIS (Koyama et al. 
2007) data from version 2.0.6.13 of the Suzaku data pro- 
cessing pipeline. Data obtained with 3x3 and 5x5 editing 
modes for each pointing were merged and standard filters 

1 We assume a flat, ACDM cosmology with Hq = 70 km s — 1 
Mpc -1 , Siju = 0.3 and fl\ = 0.7. One arc minute corresponds 
to 73 kpc at the cluster redshift (z = 0.063) in this cosmology. 
Unless otherwise noted, quoted errors correspond to 90% confi- 
dence intervals. 
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Fig. 1. Mosaic of five pointings (0.5-2 keV), with North 
at top and East at left. The image has been adaptivcly 
smoothed. Cluster spectra are taken from the solid green cir- 
cle and annular sectors. Background regions are the dashed 
blue trapezoids. Magenta circles and ellipses surround ex- 
cluded point sources. 

(requiring earth elevation greater than 5 degrees, sunlit 
earth elevation angle greater than 20 degrees, cutoff rigid- 
ity greater than 6 GeV/c, and excluding data obtained 
during and just after passages through the South Atlantic 
Anomaly) were applied. Events in detector regions nomi- 
nally illuminated by the XIS calibration sources were also 
excluded. The net good exposure times after filtering were 
10.2 ks, 23.9 ks, 24.6 ks, 27.3 ks and 38.3 ks for the central, 
near-north, near-south, far-north and far-south pointings, 
respectively. Thus the total good exposure time was 124.3 
ks. 

Smoothed images of each field were inspected by eye 
and regions centered on nine obvious point sources were 
excluded, as indicated in Figure 1. 

2.2. Image Analysis 

To obtain a properly normalized cluster surface bright- 
ness profile from our multiple pointings, we used the fol- 
lowing procedure to generate exposure maps. First, sim- 
ulated monochromatic photon lists representing a uni- 
form, 20' radius extended source were constructed, us- 
ing 10 7 photons at a number of energies for each detector 
and pointing. The energy bins were selected to repre- 
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sent roughly constant regions of the instrumental effective 
area curve. Each simulated photon list was used as input 
to the XRT ray-tracing simulator (xissim; Ishisaki et al. 
2007), which includes the effects of the optical blocking 
filter (OBF) contamination, producing an image for each 
detector, pointing, and energy combination. To increase 
signal-to-noise, these flat field images were smoothed by 
a two dimensional Gaussian kernel with 72" FWHM. The 
calibration source regions and regions of low effective area 
near the chip edges were clipped. Spectral weights were 
determined from a spectral model fit to the cosmic X- 
ray background (described in Section 2.3.2 below). These 
spectral weights were chosen under the null hypothesis 
that surface brightness fluctuations in the cluster outskirts 
are due to cosmic background variations. When applied to 
cluster emission, they will result in an overestimate of the 
cluster surface brightness, however we estimate this to be 
less than 10% (5%) in the 0.5-2 keV (2-8 keV) range for 
the cluster temperatures and abundances measured in this 
work. Errors in surface brightness between cluster regions 
of different temperatures (ranging from 2-6 keV) will be 
less than 2% due to the spectral weighting. For each chip 
and pointing, the flat field frames were combined using the 
spectral weighting to produce an exposure map. Finally, 
the individual exposure maps were mosaicked onto a com- 
mon sky coordinate (WCS) reference frame and summed. 

Images from each pointing toward the cluster were 
extracted in the 0.5 2 keV and 2 8 keV bands, and 
these were mosaicked in sky coordinates. Cutoff-rigidity- 
weighted particle background rates were computed for 
each spectral band, sensor and pointing from the archive 
of 'night-earth' data (i.e., data obtained during orbit night 
while Suzaku pointed at the earth; Tawa et al. 2008). A 
(spatially uniform) particle background component was 
then subtracted from each cluster image, and each re- 
sulting image was then divided by the appropriate ex- 
posure map. The exposure-corrected, 0.5-2 keV mosaic 
image is shown in Figure 1. For purposes of this Figure, 
the background- and exposure-corrected image has been 
smoothed with an adaptive, circular kernel. The kernel 
radius is chosen to include at least 25 counts, up to a 
maximum radius of 10 binned pixels (84"). 

2.3. Spectral Analysis and X-ray Background 

2.3.1. Analysis and Cluster Modeling 

Pulse-height spectra were accumulated for each in a se- 
ries of annular sectors 2.5-4' in width, depending on ra- 
dius. The spectral extraction regions are illustrated in 
Figure 1. Two ancillary response functions (ARF files) 
were constructed for each region and XIS detector us- 
ing the Monte Carlo FTOOL 2 xissimarf gen v2008-04- 
05 (Ishisaki et al. 2007). The first ARF was produced 
from a 20' radius uniform surface brightness source; this 
response file was used to fit uniform background emis- 
sion in the spectral analysis that follows. The second 
ARF was produced from the cluster surface brightness 
distribution results described in Section 4.1; this response 

2 http:/ /heasarc. gsfc.nasa.gov/ftools 
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was used to fit cluster emission in each region. These 
files properly account for energy-dependent vignetting and 
scattering properties of the Suzaku XRT, as well as for 
the time-dependent contamination layer on each of the 
XIS optical blocking filters. Corresponding detector re- 
distribution functions (RMF files) were constructed us- 
ing the FTOOL xisrmf gen v2007-05-14. Cutoff-rigidity- 
weighted particle-induced background spectra were deter- 
mined from night-earth data for each sensor and pointing 
using the FTOOL xisnxbgen v2008-03-08. These back- 
ground spectra were subtracted from each source spec- 
trum before fitting. 

XSPEC vl2.5.0 (Arnaud 1996) was used to fit data 
from all four XIS sensors simultaneously for all spec- 
tral fits reported here. Spectral fitting was confined to 
the 0.55-8 keV band where the response calibration is 
best. Calibration observations directly constrain the time- 
dependent, on-axis absorption by contamination of the 
XIS optical blocking filter in this spectral band (Koyama 
et al. 2007). 

The emission from each cluster region was modelled by 
simultaneously fitting the pulse-height spectra from the 
cluster region and two background regions. The sum of a 
thin thermal plasma model (mckal in XSPEC) plus a cos- 
mic background model (described in Section 2.3.2 below) 
was fit to each cluster spectrum, while only the back- 
ground model was fit to the background regions. In all 
fits the cluster heavy element abundance, relative to cos- 
mic values (Anders & Grevesse 1989), was allowed to vary. 
The best-fit redshift determined from the central 5' region 
(z = 0.0626 ± 0.0008) is consistent with the optical clus- 
ter redshift (z = 0.0625 ± 0.0003) reported by Smith et al. 
(2004). The redshift was fixed at the optical value for 
other spectral fits. The combined XIS0,2,3 (FI) spectra 
and best-fit model components for the central 5' and the 
northern extension are shown in Figure 2. 
2.3.2. Cosmic Background Model 

Spectra of the celestial X-ray background outside the 
cluster were obtained from two trapezoidal background 
regions at the northern and southern extremes of the ob- 
served field, 2.0 to 2.3 Mpc from the cluster center; see 
Figure 1. The background spectral model for each fit 
included the sum of absorbed power-law and soft ther- 
mal components, to account for extragalactic and Galactic 
backgrounds, respectively. When modelling the spectra of 
the faintest cluster regions (those with r > 17.5') we also 
included a suite of emission lines to account for possible 
contributions from geocoronal solar wind charge exchange 
(SWCX). We find that the contribution of SWCX emis- 
sion is uncertain and, though small, may be significant 
in the interpretation of the emission detected in these re- 
gions. In this section we describe the background model 
used for the brighter (r < 17.5') cluster regions, for which 
SWCX is neglected. We discuss our treatment of SWCX 
emission in the cluster outskirts in the following section. 

Best-fit spectral model parameters for the background 
regions, neglecting SWCX emission, are shown in the up- 
per half of Table 1. Consistent results (within ±6% in 
normalization) were obtained for each model component 
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Fig. 2. Spectra of the central 5' radius region and annuli toward the north from Abell 1795. The locations of the regions are shown 
in Figure 1. For each region, the black points represent the combined count rate from the three FI devices (XIS0,2,3), while the solid 
lines show contributions from individual model components: the Galactic background (blue); the extragalactic background (green); 
the cluster (red); and solar wind charge exchange (cyan). The total model is given by the black line. 



Table 1. Background Model Parameters" 



Component 



North 



South 



Joint 



without SWCX 
Power Law l\ 
Thermal kT c 
Power Law B' 
Thermal B d 
Total B d 



keV 



11.5 
0.15 
7.7 
2.9 
10.6 



+0.6 
-0.7 
+0.03 
-0.04 
+0.4 
-0.4 
+0.7 
-0.7 
+0.7 
-0.7 



10 3 +°- 6 

+0.04 
-0.02 
+0.4 
-0.4 
+0.7 
-0.6 
+0.7 
-0.7 



0.19 
7.0 
3.0 

10.0 



+0.4 
-0.4 
+0.02 
-0.03 
7 q +0.3 
'■° -0.3 
+0.5 
0.5 

10-3 ±8j 



10.9 
0.18 



3.0 



with SWCX 
Power Law l\ keV 
Thermal kT c 
Power Law B d 
Thermal B d 'f 
SWCX B d } 
Total B d 



11.3 



+0.6 

-0.7 

0.18 e 

7-6 t° .t 

2.8 tH 

1 5 + 2 ' 9 
1 -° -1.0 

11 9 +°' 7 

-0.8 



+0.7 
-0.7 

0.18 e 

+0.5 
-0.5 

2.8 ± 



10.1 



G 



i.i 

-2.2 
+ 2.3 
0.7 



1.3 

11 9 + - 7 
11 - z -0.7 



Notes: 

a uncertainties are 90% confidence for 1 interesting parameter 

b intensity, ph s _1 cm~ 2 kcV~ 1 sr _1 at 1 keV; photon index fixed at Y = 1.4 

c keV; abundance fixed at solar 

d surface brightness, 10~ 12 erg s — 1 cm -2 deg -2 , 0.5-2 keV 
e kT was fixed at 0.18 keV for the model including SWCX 
f the best-fit thermal and SWCX fluxes are anti-corrclated 
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in the two background regions. Since our spectral analy- 
sis regions are quite small, expected Poisson fluctuations 
in the extragalactic background are actually larger than 
±6%, so this level of agreement must be regarded as for- 
tuitous. To allow for these fluctuations, the normalization 
of the corresponding spectral component was allowed to 
vary separately when fitting each region of the cluster. 
This point is discussed in detail in Section 3.4 below. 

The 2-10 keV surface brightness of the power-law com- 
ponent, determined from the joint fit to both the north 
and south background regions, is B2-10 = 21.7 ± 1.5 x 
10~ 12 erg s" 1 cm -2 deg -2 , in excellent agreement with 
the value reported by DeLuca & Molcndi (2004) from 
XMM-Newton observations (22.4 ± 1.6 x 10~ 12 erg cm" 2 
s" 1 deg" 2 ). As a check on the estimated 0.5-2 keV 
background, we determined the ROSAT All-Sky Survey 
(RASS) surface brightness in an annulus centered on Abell 
1795 with inner and outer radii of 0.5 and 1 degree, respec- 
tively. We used the X-ray Background Tool 3 (Snowden 
ct al. 1997) provided by NASA's Goddard Space Flight 
Center to extract the ROSAT spectrum for this region. 
We then fit the recommended X-ray background model 
(Kuntz & Snowden 2000) , which includes two soft thermal 
components and a power law, to the ROSAT spectrum in 
the 0.1-2 keV band. This yields a surface brightness in 
the 0.5-2 keV band of Bo.5-2 = 11-1± 2 x 10" 12 erg s" 1 
cm" 2 dcg~ 2 , in excellent agreement with our (more pre- 
cise) Suzaku estimates listed in Table 1. 

Thus, even ignoring possible contributions of SWCX 
emission, our Suzaku background estimates are in excel- 
lent agreement with previous measurements of both the 
Galactic and extragalactic components. In our chosen 
spectral bands we can determine background (and fore- 
ground) flux in the vicinity of Abell 1795 with a relative 
statistical uncertainty of just over 5%. As we discuss in 
section 3.4 below, the accuracy of our background subtrac- 
tion is in fact limited by cosmic variations in the density of 
unresolved background sources rather than by statistical 
or instrumental systematic errors. 
2.3.3. Solar Wind Charge Exchange Model 

Heavy ions in the solar wind undergo charge exchange 
with neutral atoms in the Earth's gcocorona and in the 
solar magnctosphere, resulting in a field-filling emission 
line spectrum which contaminates the more distant cluster 
emission. Lines typically produced by solar wind charge 
exchange (SWCX) include those of highly-ionized C, O, 
Ne, and Mg (e.g., Snowden et al. 2004; Fujimoto et al. 
2007) . Variations in the solar wind can produce variations 
in the geocoronal SWCX flux on timescales of seconds, 
while heliosphcric SWCX contributes a stable or slowly- 
varying X-ray background (e.g., Cravens 2000). 

We have checked for SWCX variability by constructing 
light curves for each pointing, excluding the center point- 
ing where the cluster emission is expected to dominate 
any SWCX emission. Light curves were extracted for all 
four XIS detectors in the 0.4-2 keV range from the full de- 
tector field, excluding point sources. Particle background 

3 http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/xraybg/xraybg . 
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Fig. 3. (top) Combined, background-subtracted light curve 
for all detectors in the 0.4-2 keV range. The four outer point- 
ings are shown (nn: near north; fn: far north; ns: near south; 
fs: far south), (bottom) Solar wind ion flux from the ACE 
satellite, corrected for travel time to the Earth. Count rate 
variation and significant enhancements in the solar proton and 
ion flux indicate possible SWCX contamination. The green 
arrows indicate a strong count rate flare. 

light curves were constructed from the 'night-earth' data 
archive (Tawa et al. 2008) using the appropriate cut-off 
rigidity time series, and these were subtracted from the 
individual detector light curves, which were then com- 
bined for each field. The resultant complete light curve is 
shown in the top panel of Figure 3. Any residual varia- 
tions in count rate are a result of variations in the emission 
from the Galactic foreground, the unresolved extragalac- 
tic background, the cluster, and/or SWCX. The extrac- 
tion regions for corresponding pointings are identical in 
sky coordinates, so any differences between far north and 
far south ("fn" and "fs") or between near north and near 
south ( "nn" and "ns" ) must be due to variations (tempo- 
ral or spatial) in the surface brightness of these emission 
components. 

There is clear variability in the light curve. Most 
notably, the far north count rate shows a strong spike 
(marked by the green arrow). Count rates in the other 
three pointings are also inconsistent with the null hy- 
pothesis of constancy (p- values < 0.1), although they 
lack strong variation. Analysis of public data from the 
Advanced Composition Explorer (ACE) satellite (Stone 
et al. 1998) reveals an enhanced flux of protons and heavy 
ions during much of the observation (see Figure 3, bot- 
tom panel). Although there are gaps in the ACE data, 
the proton, C +6 , and +7 flux during the far north and 
near north pointings (near MJD 53716) arc similar to lev- 
els observed during times of SWCX contamination in a 
Suzaku observation of the North Ecliptic Pole (Fujimoto 
,gt al. 2007). The far north count rate spike occurs near a 
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Table 2. Modeled SWCX Emission Lines 



transition 


energy 


photon flux 




(keV) 


(ph s _1 cm~ 2 sr _1 ) 






JNortn 


boutn 


vn triplet 


0.574 


o no +4.U7 
3 - 03 -1.56 


< 6.00 


vili Lya 


0.654 


< 2.25 


< 1.77 


z' \ - rTTX T /D 

U VIII Lyp 


0. 114, 


< 0.71 


< 0.61 


vin Ly7 


0.817 


< 0.63 


0.49 


Ne ix triplet 


0.923 


< 0.47 


4fi + - 38 
u -^° -0.35 


Ne x Lya 


1.022 


< 0.31 


41 +°' 28 


Mg xi triplet 


1.345 


< 0.11 


< 0.23 


Mg xn Lya 


1.472 


n 4s +°- 23 

U -^ J -0.22 


< 0.16 


Mg xn Ly/3 


1.745 


< 0.08 


< 0.11 



peak in both the proton and heavy ion flux. There is an 
additional spike in the heavy ion flux near the end of the 
near north observation. The solar particle flux near the 
end of the observations (after MJD 53716.5) is consistent 
with that observed during quiescence in other Suzaku ob- 
servations (Fujimoto et al. 2007; Miller et al. 2008). We 
conclude that our observations are likely contaminated by 
variable SWCX emission. 

We estimate the level and variability of gcocoronal 
SWCX contamination by comparing the emission from 
the (presumed cluster-free) background regions. As in 
Section 2.3.2 above, the non-SWCX background is mod- 
eled as the sum of a soft thermal (Galactic) background 
component and a power-law (extragalactic) component. 
When fitting the spectra of the relatively small back- 
ground regions alone, we found strong degeneracies be- 
tween the SWCX and thermal model parameters, so we 
fixed the temperature of the thermal model at kT = 0.18 
keV, the value determined from the joint SWCX-free 
model. The photon index of the power-law component 
was fixed at T = 1.4. The normalizations of both back- 
ground components were allowed to vary independently 
in each region (see Section 3.4). 

To model the SWCX in the background regions, we in- 
clude a suite of unresolved Gaussian emission lines typ- 
ically observed in SWCX emission between 0.5-2 keV 
(Snowden et al. 2004; Fujimoto et al. 2007), listed in Table 
2. The line energies are fixed to their theoretical values, 
using the energy of the dominant forbidden transition for 
helium-like lines of O vil, Ne IX, and Mg xi (Koutroumpa 
et al. 2007). The line fluxes arc allowed to vary indepen- 
dently. 

We use Markov chain Monte Carlo (MCMC) techniques 
to constrain the relatively large number of model param- 
eters and their uncertainties. The MCMC tools in XSPEC 
vl2.5.0 (Arnaud 1996) were used to produce 8 chains of 
7000 steps, of which the first 2000 were excluded as "burn- 
in" . The input covariancc matrix was scaled to get suf- 
ficient mixing (i.e., so that the fraction of repeated loca- 
tions in parameter space was ~ 0.8). These 40,000 chain 
steps were combined to produce probability distributions 
for each of the 22 free parameters, and from these the best 
fit value and 90% confidence intervals were estimated. 



Both the far north and far south background regions 
show an excess of line emission above the Galactic and ex- 
tragalactic components. The spectra in Figure 4 suggest 
enhancement of the O vn triplet near 0.56 keV (evident 
in XIS1) and a ~ 3-cr feature near 1.5 keV (evident in 
XIS0,2,3) in the far north spectrum. The latter is consis- 
tent with Mg xn Lya. The far south spectra show much 
weaker upper limits to the enhancement near 0.56 keV (in 
XIS1) and 1.5 keV (in XIS0,2,3), but they contain excess 
emission in the 0.9-1.0 keV region (evident in all detec- 
tors) that is consistent with Ne IX and Ne x. The line 
strengths and upper limits arc shown in Table 2. The 
Mg xn line is similar in energy to a strong instrumental 
Al Ka line, however the Al Ka line flux would have to be 
enhanced by more than a factor of 3 during the full far 
north exposure to produce such strong residual emission. 
Background variations of this sort have not been reported 
in previous Suzaku observations, nor in the night-earth 
background analysis of Tawa et al. (2008). The presence 
of a Mg xn SWCX line without additional SWCX lines 
such as O vin, Ne ix and Ne X seems puzzling, although 
since detection of Mg xn with XMM-Newton is compli- 
cated by a much stronger and more variable instrumental 
Al Ka line, there are few data from other sources. Mg XII 
Ly7 has likely been observed in another XMM-Newton 
SWCX event (Carter & Sembay 2008). Snowden et al. 
(2004) claim a Mg xi SWCX detection at 1.34 keV. The 
upper limits on the photon fluxes for other expected lines 
are low compared to other SWCX detections, which sug- 
gests that our observations are less affected by geocoronal 
SWCX. 

The contributions of SWCX and other model compo- 
nents to the total surface brightness in the background 
regions are listed in the lower half of Table 1. The Table 
shows that SWCX emission is detected at better than 
90% confidence during our observations of both far north 
and far south fields, and that the SWCX emission spec- 
trum differs significantly between the two observations. 
However, in each field SWCX contributes less than 15% 
of the total, soft-band cosmic background brightness, and 
the difference between SWCX contributions to the two 
fields is less than 5% of the total cosmic background 
brightness. 

Our spectral analysis of the outermost regions of the 
cluster (r > 17.5') follows that used for the interior of 
the cluster, except that we include the SWCX model de- 
scribed above in the background model (see Figure 2). In 
particular, the soft thermal (Galactic) background com- 
ponent was fit with a single free temperature and nor- 
malization for all regions (17.5-26' north and south, and 
the north and south background regions). The power-law 
(extragalactic) background component was fit with a fixed 
photon index T = 1 .4 and normalizations that were allowed 
to vary independently in each region (see Section 3.4). A 
mekal component, with temperature, abundance, and nor- 
malization free to vary, was used to model cluster emis- 
sion. 



No. 



Suzaku Observations of Abell 1795 



XISl fn 



M i l fa - .-^-d^dv; 



,± +r- 



XIS1 fs 



I ^ r ■!, 1 ■ i i I + 



fe^tj^tt vt^ ^* + ,* + + >+ 4^-i*^ 



XISl fn 
XISl fs 



0.5 



1 1.5 
energy (keV) 



XIS0,3,3 fn 



XISO.2,3 fs 



XIS0,2,3 fn 
XISO.2,3 fs 



0.5 



1 1.5 
energy (keV) 



Fig. 4. Spectra of the two background regions from the BI detector XISl (left) and the combined FI detectors XIS0,2,3 (right). The 
labels "fn" and "fs" refer to the far north and far south pointings, respectively. In the top panels, the red lines indicate the best-fit 
diffuse Galactic + extragalactic background, and the cyan lines indicate 90% upper limits to SWCX emission modeled by Gaussian 
lines (see Table 2). The residuals are calculated from the combined background + SWCX model. The model was fit simultaneously 
to spectra from the four detectors; the FI spectra were only combined here to improve the displayed S/N. The bottom panels directly 
compare the counts spectra from the far north (black) and far south (red) background regions. Differences can be seen near the 
O vii (0.56 keV), Ne ix and Ne x (1 keV), and Mg xn (1.5 keV) lines. 



3. Systematic Errors 

We wish to measure cluster emission at very low sur- 
face brightness levels and it is important to understand 
the influence of various instrumental and other systmatic 
effects that may limit our ability to do so. In this sec- 
tion we consider a number of potential systematic errors 
in more detail, including i) the impact of our choice of 
passband on temperature measurement accuracy; ii) the 
accuracy of our knowledge of spatial response variations 
over the XIS field of view; iii) the impact of scattered light, 
both from the bright cluster core and from point sources 
in the field of view, on the accuracy of our spectral mod- 
elling; and iv) the effects of counting statistics on the flux 
of unresolved sources. 

3.1. Spectral passband, Galactic absorption and cluster 
temperature 

A drawback of our decision to ignore data at energies 
below 0.55 keV in our spectral fits is that it is difficult 
to constrain the Galactic absorption directly from the X- 
ray data unless the signal-to-noise ratio is quite high. On 
the other hand, we have found that the best-fit model 



temperature and Galactic column do exhibit a mild (anti) 
correlation. Thus, any systematic errors in determination 
of the Galactic column (such as might arise from spatially 
varying errors in the contamination correction) can induce 
errors in the measured temperature. We have therefore 
elected to fix the model Galactic column density at the 
value obtained by fitting to the high-signal-to-noise ratio 
spectrum of the integrated cluster emission at r < 5.0'. 
Fitting results are shown in Table 3. The best-fit value 
of uh is in excellent agreement with the nominal Galactic 
column (riff = 1.2 x 10 20 enr 2 ) estimated for this field by 
NASA's HEARSARC tool nH. 4 

3.2. Spatial response variations and reproducibility of 
spectral modelling 

As a check on the reliability of the Suzaku calibration, 
we compared fitting results for regions of the cluster ob- 
served with multiple pointings. Two such 'overlap' regions 
provide sufficient flux for a meaningful comparison. Each 

4 http: //heasarc/gsf c .nasa. gov/cgi-bin/Tools/w3nh. The 

column densities estimated from the Leiden/Argentine/Bonn 
(Kalberla et al. 2005) and Dickey & Lockman (1990) maps agree 
within 1%. 
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Table 3. Best- fit Model Parameters for r < 5.0' 



Flux a 


kT b Abundance 




X 2 /D0F 


7 46 +(U)4 


5.26 ±0.06 0.360 ±0.018 


1.2 ±0.3 


4709.6/4481 


Notes: 








a 10~ n erg 


s- 1 cm" 2 , 0.55-8 keV; b keV; c 


10 20 cm' 2 





Near-north Pointing (+) 
Center Pointing (x) 




0.01 



0.02 0.03 
n„(10 22 cm- 2 ) 



0.04 



0.05 



Fig. 5. Contours of x 2 over the rijj-kT plane derived from 
two distinct observations of region of Abell 1795. The re- 
gion in question is an annular sector 5'— 7.5' north of the clus- 
ter. Contours at 68% (black), 90% (red) and 95% (green) 
confidence from the 'near-north' (solid contours) and 'center' 
(dashed contours) pointings are indicated. The cross and x 
symbols indicate the best fit values for these two regions, re- 
spectively. The good agreement of these measurements from 
opposite sides of the XIS field of view illustrates the accuracy 
with which spatial response variations have been calibrated. 

region is an annular sector with inner and outer radii of 
5.0' and 7.5', respectively. The northern overlap region 
was observed with central and 'near-north' pointings; the 
southern overlap region was observed with the central and 
'near-south' pointings. Thus each overlap region was ob- 
served on opposite sides of the XIS field of view. In order 
to gauge errors in the calibration of the soft response, 
we allowed the Galactic column to vary when fitting indi- 
vidual spectra. For each overlap region, the best-fit model 
parameters from the two pointings are consistent with one 
another within 90% confidence statistical errors (SkT = 0.3 
keV). The level of agreement, as well as the magnitude of 
calibration errors, is illustrated in Figure 5, which shows 
contours of x 2 m the- nn-kT plane for the northern over- 
lap region. Figure 5 suggests residual systematic errors in 
the determination of the absorbing column may be of or- 
der 5nn ~ 1 x 10 20 cm" 2 . As noted above, in subsequent 
analysis we fix uh = 1-2 x 10 20 cm" 2 and include a sys- 
tematic temperature error of 5kT = 0.15 keV to account 
for this effect. 

3.3. Scattered X-ray flux 

We estimated the effect of X-rays scattered by the 
XRT from the bright, cool core of Abell 1795 to the 



cluster outskirts with the aid of a number of simula- 
tions. For this purpose we used the FTOOLS xissim and 
xissimarf gen (Ishisaki et al. 2007). For fields centered 
20' from a bright source (comparable to our far north and 
far south fields), actual stray intensity levels are expected 
to be less than ~ 2 times the values predicted by these 
tools (Serlemitsos et al. 2007). 

To gauge the effect of scattered X-rays from the bright 
cluster core on the surface brightness profiles at large ra- 
dius, we simulated observations of a point source located 
at the cluster center with a spectrum and total flux ob- 
served within 5' of the center (see Table 3). Since the 
actual cluster emission is more diffuse than assumed in 
this model, our simulation overestimates the scattered flux 
from the (unresolved) bright cool core of the cluster. We 
simulated the full mosaic observation of five pointings for 
each sensor, and normalized the simulated observations 
using the exposure maps discussed in Section 2.2 above. 
We then extracted surface brightness profiles for this point 
source using the same methods we used to extract the 
cluster profiles (as discussed in Section 4.1). In particu- 
lar, all regions excluded from the cluster profiles because 
they contained cosmic sources were also excluded in the 
simulations. For clarity, we did not add a cosmic X-ray 
background signal to the simulated point source profile. 

The resulting point-source surface brightness profile in 
the 0.5-2 keV band is shown by the points in Figure 6. 
The Figure shows that the expected scattered flux at 
r > 10' exhibits a distinct oscillation at relatively low am- 
plitude. The oscillation is presumably due to the vari- 
ations with field angle of the intensity of so-called 'sec- 
ondary' (single bounce) and 'backside' (three-bounce) re- 
flections (Serlemitsos et al. 2007). 

Comparison of Figure 6 with the estimated celestial 
background brightness listed in Table 1 shows that the 
surface brightness due to scattered flux in this radial range 
is always less than 7% of the cosmic background surface 
brightness we observe at the largest radii (r > 30') . The ac- 
tual scattered intensity from the cluster will be below that 
shown in Figure 6 at these radii because convolution of the 
point response function (PRF) with the extended cluster 
flux distribution will reduce the amplitude of the oscillat- 
ing wings. Serlemitsos et al. (2007) note that, given align- 
ment tolerances of the XRT reflectors and pre-collimator 
blades, the actual scattered flux at r = 20' could in princi- 
ple exceed the ray-tracing simulations (which assume per- 
fect alignment) by as much as a factor of 2. Allowing for 
this uncertainty, we expect that scattered cluster flux will 
contribute at most 10% of the total observed brightness 
at r > 10'. 

We confirmed the effect of scattered cluster flux on 
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Fig. 6. Simulated surface brightness (points) in the 0.5—2 
keV band for a point source located at the cluster center with 
flux equal to that from within r < 5' of the cluster center. 
Surface brightness northward (black circles) and southward 
(red squares) from the source location are shown. No cosmic 
background has been added. 

our spectral analysis in the outermost regions of our field 
(17.5' < r < 32') by computing the effective area (as a 
function of energy) for detection of photons there origi- 
nating from the central 5' of the cluster. Here we used 
the measured cluster surface brightness distribution as in- 
put to the effective area calculation, rather than a sim- 
ple point source. Note that this central region produces 
more than 80% of the total cluster flux reported for A1795 
in the Northern ROSAT All-Sky Survey (NORAS) cata- 
log (Bohringcr ct al. 2000). Using these effective area func- 
tions, and the cluster spectral model described in Table 3, 
we employed XSPEC to simulate the flux and spectrum of 
X-rays scattered from the cluster center into each cluster 
(and background) region at r > 17.5'. The resulting fluxes 
are listed in the first two columns of Table 4. 

We conclude from Table 4 that the scattered intensity 
from the cluster core contributes ~ 3-5% of the total ob- 
served surface brightness in the regions for which we per- 
form spectroscopy. Both magnitude and variation with 
radius are in agreement with the foregoing analysis of the 
PRF. We reiterate that uncertainties in collimator align- 
ment render these estimates uncertain by as much as a 
factor of 2. The estimated scattered fluxes in the north- 
ern and southern regions differ somewhat, particularly in 
the radial range 17.5-22'. The differences are presumably 
due to the exclusion of point sources in the north (see 
Figure 1), which changes the effective radii of the annular 
bins there, coupled with the rapidly varying scattered flux 
in this region (see Figure 6). In any event, the scattered 
flux in all regions listed in Table 4 is comparable to or less 
than the statistical errors in the estimated background 
(see Table 1). We note that, because of the oscillating 
PRF wings (see Figure 6), the simulations imply that spa- 
tial non-uniformity of the scattered flux will cause us to 
over-estimate the background at r > 27.5', and thus to 



(slightly) underestimate the cluster flux in the outer-most 
regions. 

Finally, we investigated the effect of scattered flux from 
bright point sources on our cluster spectrophotometry. 
The brightest such source is in the far-north field, with 
Suzaku position 5 a = 13h49m03.0s,<5 = +26°58'46.1", and 
we estimate its flux to be 4.0 ± 1.9 x 10 -13 erg s _1 enr 2 
(0.5-2 keV) within an aperture of 2.5' radius. We have ex- 
cluded a 3.5' diameter region around this source from the 
cluster analysis. We estimated the flux contributed by this 
source to each of our spectroscopic regions following meth- 
ods described above for the cluster flux. Results are listed 
in Table 4. The wings of the image of this bright point 
source contribute less than 5% of the observed total flux 
to each of the two outermost spectral analysis regions (ex- 
tending from 17.5' to 26' from the cluster). Thus in these 
two regions, the bright point source contribution is com- 
parable to that from scattered cluster flux. In all other 
regions, the bright point source contribution is smaller. 
We note that point source contributions should be known 
more accurately than those from the cluster core because 
the the ray-tracing simulations are better calibrated at 
r ~ 3'-10' than at r > 20' (Serlemitsos et al. 2007). 

3.4- Cosmic background variations 

Wide and deep X-ray surveys conducted over the past 
~ 15 years have revealed a great deal about the distribu- 
tion in flux and space of the sources that constitute the 
extragalactic X-ray background (Brandt & Hasinger 2005 
and references therein) . We have used this information to 
estimate the expected variation of the extragalactic back- 
ground not resolved by Suzaku. On scales larger than 
the Suzaku PRF (~ 2'), the correlation of source angu- 
lar positions is weak (Yang et al. 2003), and for purposes 
of estimating the expected variance in our background, 
we neglect it. Formally the expected background sur- 
face brightness B due to sources with flux S < S exc i is 
B = ^" cl ^SdS, where is the differential number of 
sources per unit solid angle and flux. The expected vari- 
ance of the background brightness, when measured 

over a solid angle f2, is a% = {J^ sxcl ^S 2 dS)/n. We 
adopt a broken power-law model for ^g using parameters 
derived by Moretti ct al. (2003). 

It may be instructive to note that if the integral number 
counts N(S > S exc i) varied simply as N(S > S exc i) oc S~" cl , 
then the relative fluctuation in the the unresolved back- 
ground within a solid angle Q would satisfy ctb/B = 
K(a)/{N(S > Sexd)ty 1/2 , where the factor K(a) varies 
between ~ 1/2 and 3/2 for appropriate values of a. This 
relationship sets the flux limit S exc i to which point sources 
must be detected for accurate background estimation in a 
field of solid angle CI. 

Our Suzaku data alone allow us to detect point sources 
to a limiting flux S > S exc i = 105*14, where Si 4 is the flux in 
units of 10~ 14 erg s _1 cm~ 2 , in cither the soft (0.5-2 keV) 
or hard (2-8 keV) bands. For a significant fraction of our 

5 This source is coincident, within expected Suzaku position error 
of up to 20 arcsec, with 1RXS J134903.6+265845. 
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Table 4. Estimated scattered flux 



From Cluster Core From N. Point Source Total 



Region a 


r>b 

sc. 0.5- 


-2 sc. 2 — 8 


S C , U . O — z 


Di) 

sc,2 — 8 


R& 

SC,U.t> — ^ 


sc. 2 — 8 


17.5-22 (N) 


0.36 


0.25 


0.38 


0.02 


0.74 


0.27 


17.5-22 (S) 


0.21 


0.16 






0.21 


0.16 


22-26 (N) 


0.37 


0.33 


0.45 


0.03 


0.82 


0.36 


22-26 (S) 


0.37 


0.36 






0.37 


0.36 


27.5-31.7 (N) 


0.47 


0.35 






0.47 


0.35 


27.5-31.7 (S) 


0.45 


0.34 






0.45 


0.34 


Notes: 

a Minimum and 
b lCT 12 erg s- 1 


maximum 
cm -2 deg" 


radii, arcmin, north (N) or 
" 2 , 0.5-2 and 2-8 keV 


south (S) of cluster center. 





field, we can also use archival XMM-Newton observations 
to constrain the number of unresolved sources contribut- 
ing to the Suzaku background at fluxes as low as S exc i = 
IS14 in the soft band, and S exc i = 1.3 Si 4 in the hard band 
(see below) . Using our adopted source-count flux relation, 
we find the expected RMS fluctuations in background due 
to unresolved sources are ctb = 3.9(0.94) x 10 _12 £7q q{ 2 erg 
s _1 cm -2 deg -2 for S exc i = 10(1) Su in the soft band. The 
corresponding fluctuation amplitudes in the hard band 
are <jb = 4.5(2.5) x 10~ 12 fi^ g{ 2 erg s -1 cm -2 deg~ 2 for 
Sexci = 10(1.3) S14. Here fio.oi is the solid angle of the 
measurement region in units of 10 -2 deg 2 . 

The field coverage of the available XMM-Newton ob- 
servations is shown in Figure 7. We searched the Second 
XMM-Newton Serendipitous Sky Survey (XMM2P) cata- 
log available at NASA's High Energy Astrophysics Science 
Archive Research Center 6 for sources in our Suzaku field 
at distance 10' < r < 26' from the cluster center. We 
found that the effective (EPIC PN) exposure in this re- 
gion ranges from 10-50 kscc. Our adopted XMM-Newton 
flux thresholds (S exc i = IS14 and S exc i = 1.3 Su in the soft 
and hard bands, respectively) exceed the nominal XMM- 
Newton 5cr detection thresholds (Watson et al. 2001) for 
the minimum effective exposure available in our field. We 
summed the flux of sources brighter than these thresh- 
olds in each of our spectral analysis regions. Since the 
XMM2P catalog does not tabulate source flux in the 2-8 
keV band, we estimated the 2-8 keV flux by scaling the 
tabulated XMM-Newton 2-12 keV flux by a factor ap- 
propriate for power-law photon number distribution with 
index 1.4 (that is, we assumed S2-8 = 5*2-12 x 0.67.) The 
contributions of sources resolved by XMM-Newton to the 
Suzaku background, as well as the estimated fluctuations 
in the background due to fainter sources, are listed in 
Table 5. 

In the outskirts of the cluster (r > 10'), the expected 
fluctuations in the unresolved background arc ~ 10 -12 
erg s -1 cm -2 deg -2 (0.5-2 keV), or about 10% of the 
total (Galactic plus extragalactic) background. In our 
two background regions at r > 27.5', which lack complete 




26:00:00.0 I -I 

30.0 13:49:00.0 30.0 48:00.0 

Rjght ascension 

Fig. 7. XMM-Newton image of Abell 1795 field with Suzaku 
spectral analysis regions superposed in green. Where avail- 
able, the XMM-Newton data constrain the background un- 
certainties due to point sources unresolved by Suzaku. As in 
Figure 1, background regions are the dashed blue trapezoids 
and magenta circles and ellipses show point sources excluded 
from the Suzaku data. 
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Table 5. Estimated Background Contributions from Point Sources 



Radius" 


0.01 


CC 
exel 


0.5- 


-2 keV 

U Tl V C S 


u ri v c s 


ce 

exel 


2-1 

B L 


3 keV 

UTIVC-S 


UTIVGS 


North: 




















5-7.5 


0.60 


10 


s 


7.8 


4.9 


10 


s 


16.4 


5.8 


7.5-10 


0.85 


10 


s 


7.8 


4.2 


10 


s 


16.4 


4.9 


10-12.5 


0.93 


1 


2.2 


3.1 


0.96 


1.3 


2.7 


12.0 


2.6 


12.5-15 


0.98 


1 


3.9 


3.1 


0.94 


1.3 


3.8 


12.0 


2.6 


15-17.5 


0.91 


1 


1.3 


3.1 


0.98 


1.3 


0. 


12.0 


2.6 


17.5- 26 


1.58 


1 





3.1 


0.74 


1.3 


7.4 


12.0 


0.8 


27.5-31.7 


1.54 


10 


s 


7.8 


3.1 


10 


s 


16.4 


3.6 


South: 




















5-7.5 


0.57 


10 


s 


7.8 


5.1 


10 


s 


16.4 


6.0 


7.5-10 


0.79 


10 


9 


7.8 


4.3 


10 


9 


16.4 


5.1 


10-12.5 


0.93 


1 





3.1 


0.96 


1.3 





12.0 


2.6 


12.5-15 


0.83 


1 


1.3 


3.1 


1.0 


1.3 


4.6 


12.0 


2.7 


15-17.5 


1.22 


1 


2.0 


3.1 


0.84 


1.3 


3.8 


12.0 


2.3 


17.5-26 


3.11 


1 





3.1 


0.54 


1.3 


1.1 


12.0 


1.4 


27.5-31.7 


1.25 


10 


s 


7.8 


3.5 


10 


ff 


16.4 


2.2 



Notes: 



a Radial range of region relative to cluster center, arcmin; b Region solid angle, 10~ 2 deg 2 ; c Limiting flux 
for resolved sources, in units 10 — 14 ergs s _1 cm~ 2 , 0.5—2 keV; d 10 — 12 ergs s~ 1 cm~ 2 dcg~ 2 , 0.5-2 keV; e 
Limiting flux for resolved sources, in units 10~ 14 ergs s — 1 cm~ 2 , 2-8 keV; * 10 — 12 ergs s — 1 cm~ 2 deg -2 , 
2—8 keV; 9 These regions lack useful XMM-Newton coverage. 



XMM-Newton coverage, the expected fluctuations are sig- 
nificantly larger, ranging from 3-3.5 x 10 -12 erg s _1 cm" 2 
dcg^ 2 . We note also that the expected extragalactic back- 
ground brightness for these regions is in excellent agree- 
ment with the measured values (see Table 1). 

Finally, we stress that the expectation value of the ex- 
tragalactic background brightness varies from region to re- 
gion because the extragalactic source population (as well 
as our knowledge of it) varies (see Table 5). Thus, for ex- 
ample, it is appropriate when modelling the background 
to allow the extragalactic component to vary in brightness 
from region to region, and we do so in our fits. 

4. Results and Discussion 

4-.1. Surface Brightness 

Radial surface brightness profiles were extracted from 
exposure-corrected images in the 0.5-2 keV and 2-8 
keV spectral bands using the Chandra X-ray Center's 
CIAO/Shcrpa tools. Point sources indicated in Figure 1 
were excluded. The center of the profile was determined 
from a fit of a circular beta-model (S(r) ~ (l + j; 2 ) - ' 3 ' 3- ^' 
with x = — ) to the (two-dimensional) 0.5-2 keV image. 
A one-dimensional model consisting of a /3-modcl plus a 
constant background was then fit to each of the radial 
profiles. The profiles and the best-fitting models are pre- 
sented in Figure 8. The count-rate to flux conversions for 
these profiles were derived from the SWCX-free spectral 
model of the emission from the background regions (see 
the upper panel of Table 1). As explained in Section 2.2, 
this leads to a cluster flux overestimate of less than 10% 
(5%) for the 0.5-2 keV (2-8 keV) band. It is appropriate 



for our null hypothesis that surface brightness variations 
in the A1795 outskirts are due to cosmic background fluc- 
tuations. 

The northern and southern surface brightness profiles 
shown in Figure 8 agree remarkably well with one another 
(and with a simple /J-model with (3 = 0.64± 0.01) for r < 1.0 
Mpc (14') «r 2 5oo- As expected, the cluster is symmetrical 
and evidently relaxed in this region. At larger radii, the 
profiles diverge, with the surface brightness apparently 
falling more rapidly than the /3-model in the south, but 
actually rising in the north to an apparent peak near r = 
1.9 Mpc (26') = r 20 o in the 0.5-2 keV band. The surface 
brightness in this peak is about twice that in the nominal 
background regions, which are at r > 2.0 Mpc (r > 27.5'). 
The southern profile actually drops ~ 20% below the level 
in the nominal background regions. 

The surface brightness profiles at radii r > 10' are re- 
plotted in Figure 9, along with estimates of the cosmic 
X-ray background as a function of radius. To minimize 
the Poisson fluctuations in the extragalactic background 
we have re-binned the data at r > 18'. The background 
estimates, shown as dashed histograms in Figure 9, in- 
clude three components. We assume a spatially uni- 
form Galactic contribution in the 0.5-2 keV band, with 
flux equal to the mean thermal flux obtained from the 
spectral fitting of the background regions (see Table 1): 
B G ai =3.0 ±0.3 x 10~ 12 erg s" 1 cm" 2 deg" 2 {la errors). 
This estimate assumes no SWCX emission; we consider 
the effects of SWCX below. We adopt as the extragalactic 
background the sum of the total flux of any point sources 
detected by XMM-Newton plus an estimate of the unre- 
solved flux derived from our adopted model of number 
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Fig. 8. Surface brightness profiles northward (black circles) and southward (red squares) from the cluster center in the 0.5-2 keV 
(upper panel) and 2-8 keV (lower panel) bands. The solid curves show the best-fit /3-models to all the points (northward and 
southward) in each panel. 
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Fig. 9. Surface brightness profiles at large radius in the 0.5-2 keV (top) and 2—8 keV (bottom) bands. Profiles running north- 
ward (black circles) and southward (red squares) from the cluster center, with estimated corresponding background levels (dashed 
histograms; black to north, red to south) are shown. The solid curves are the best- fit /3-models assuming a spatially constant back- 
ground. The cluster surface brightness data have been rebinned to match the spectral extraction regions; the models are identical to 
those shown in Figure 8. The background estimates account for point sources detected by XMM-Newton but unresolved by Suzaku, 
and assume no SWCX emission (see text). 
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counts (Moretti et al. 2003). The latter two components 
are listed for each spectral analysis region in Table 5. 
Error bars plotted on the background histograms include 
the expected RMS Poisson fluctuations, also as listed in 
Table 5. The global best-fit (3- models (assuming a spa- 
tially uniform background) are shown as solid curves in 
Figure 9. 

Figure 9 shows that the background model agrees well 
with the measured surface brightness in the background 
regions (at r > 26') , which suggests that our extragalactic 
background model is reasonable. In the region 17.5' < 
r < 26' (1.3 < r < 1.9 Mpc; rsoo < r < f2oo)- the observed 
soft-band emission significantly exceeds the background 
in both the north and the south. The formal net soft- 
band cluster surface brightnesses in this radial range are 
B c i = 8.3±0.8 (2.2±0.6) x 10" 12 erg s" 1 cm" 2 deg" 2 in the 
north (south), quoting la errors. Allowing for scattered 
flux contributions listed in Table 4, the north detection 
is significant at a level greater than 8<r; the net cluster 
surface brightness in the south falls below 3er. 

If we allow for possible SWCX contributions, the esti- 
mated Galactic plus geocoronal background is somewhat 
higher, but we still find significant cluster emission in the 
north. In this case, the estimated Galactic plus geocoro- 
nal background rises to B Ga i +SWC x = 4.3 ± 0.5 (4.1 ± 
0.5) x 10~ 12 erg s" 1 cm" 2 deg -2 in the north (south), 
with 1 a errors. The corresponding results for cluster flux 
are B cl = 7.0 ±0.8 (1.1 ±0.6) x 10" 12 erg s" 1 cm" 2 deg" 2 
in the north (south), with la errors. In the north, the 
signal exceeds the background by at least 6.8 a, allow- 
ing for the scattered light contributions listed in Table 4. 
The 3<7 upper limit on the cluster emission in the south 
is B c i < 1.8 x 10" 12 erg s _1 cm -2 dcg~ 2 . 

We also note that formally the measured net cluster 
surface brightness in the north exceeds that in the south 
by a factor of 6.4 ±3.6 (la error); the absolute bright- 
ness difference of 5.9 ± 1.1 x 10~ 12 erg s" 1 cm" 2 deg~ 2 is 
significant at more than 5.3 a. We conclude that Suzaku 
has detected soft-band cluster emission with high confi- 
dence in a region with rsoo < r < r2oo north of the cluster 
center, but not in the corresponding region in the south. 
In this radial range, the cluster is significantly brighter in 
the north than in the south. 

The lower panel of Figure 9 shows that Suzaku detects 
no hard-band cluster emission at r > 17.5' (1.3 Mpc) in ei- 
ther north or south. Indeed, the Suzaku surface brightness 
measurements are remarkably consistent with the back- 
ground estimates derived from the XMM-Ncwton point- 
source catalog and the assumed number-counts model for 
17.5' < r < 26'. In the hard band, the significant north- 
south difference in measured Suzaku surface brightness 
is attributable entirely to differences in the background 
source populations in these two regions. 

Finally, comparison of the top and bottom panels of 
Figure 9 shows that our background model implies spatial 
variations in the spectral shape of the background. For 
example, the background is marginally harder in the north 
than in the south in the radial range 17.5' < r < 26', and in 
the north the background is harder at 17.5' < r < 26' than 



Table 6. Power-law fits to surface brightness profiles 



Direction 


Index" 




0.5-2 kcV 


2-8 kcV 


r > 7.5' 


N 
S 

N+S 


9 fiq+0.32 
z - uo -0.31 
o 4 q+0.29 

2 98 +0 ' 21 


o 7 q+U.44 
°- ' ^-0.38 

4 06+ 45 

q oc+0.31 
J -c.°-0.29 


7.5' <r< 17.5' 


N 
S 

N+S 


o ln +0.36 
°- lu -0.36 
q qn+0.33 
°- ou -0.32 

o 2 n+ - 24 

J - zu -0.24 


q em +0.51 
°- oi -0.47 

q fi»+ - 48 
3-68±«;It 



Notes: 



a Index a in S(r) oc r a , with 90% errors. 



at r > 26'. We allow for these variations in the spectral 
fitting used to extract the temperature profiles discussed 
below. 

Simulations of cluster X-ray emission in the vicinity of 
the virial radius (Roncarelli et al. 2006) suggest that the 
surface brightness profile follows a broken power law with 
the break at r m r2oo- Simple power-law fits to the A1795 
surface brightness profiles at r < r2oo are roughly consis- 
tent with these simulations. Best-fit values for the power- 
law index a, where S(r) oc r -Q , are listed in Table 6. The 
two different radial ranges shown correspond to r > 0.3r2oo 
(r > 7.5') and 0.3r 200 < r < 0.7r 200 (7.5' < r < 17.5'). In 
the soft band, the northern and southern profile slopes 
are identical, within errors, over the smaller radial range, 
but the apparent excess emission at r > 0.7r2oo (r > 17.5') 
flattens the profile there. In all cases the hard-band pro- 
file slopes are steeper than the corresponding soft-band 
slopes. 

The soft-band power-law index for the southern pro- 
file is quite similar to that reported by Neumann (2005) 
from stacked ROSAT profiles for a number of hot clus- 
ters. For example, Neumann finds a = 3.79^0 37 f° r 
0-3r2oo <r < 1.2r2oo for a set of clusters similar to (and in- 
cluding) A1795. On the other hand, the soft-band Suzaku 
profile north of the cluster center is somewhat shallower 
than the stacked ROSAT profile. 

We have fit the rolling power-law function used by 
Roncarelli et al. (2006) to our data. This function is 
S(x) oc x' 1 ^ with x = r/r 2 oa and y(x) = —(b max x + 
b m i n ) / (1 + x) . At small r, S(x) ^ x ~ bmin , while at large r, 
S(x) — > x~ bmal , with the transition occurring around x~ 1. 
Our data are confined to i £ 1. so they do not constrain 
b m ax well. Therefore in our fits we fix b max = 5.3 (7.5) 
for the soft (hard) band, values typical of the simulated 
clusters (Roncarelli et al. 2006). With these constraints 
we find b m in = 2.6-3.5, consistent with the simulations, 
for both hard-band profiles and for the soft-band profile 
in the south. Again because of the excess emission at 
r > 17.5' in the north, however, we find b m i n ~ 2.2 for the 
soft-band profile in this direction; this is flatter than the 
simulated clusters. 

We note that Roncarelli et al. (2006) filter the emission 
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o Suzaku N (90% errors) 
□ Suzaku S (90% errors) 
- Chandra (67% errors) 
XMM (90% errors) 




5 10 15 20 25 30 
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Fig. 10. Projected temperature profiles running northward 
(black circles) and southward (red squares) from the clus- 
ter center. Data from Chandra (green triangles; Vikhlinin ct 
al. 2006) and XMM-Newton (blue diamonds; Snowdcn ct al. 
2008; Snowdcn, private communication) are also shown. The 
radius values for the Suzaku data reflect the average origi- 
nating radius of photons detected in each extraction region, 
as described in Section 4.2. Likewise, the horizontal Suzaku 
error bars indicate 90% confidence for originating location of 
detected photons. The error bars for other missions indicate 
the radial range of the bin. 

from their simulated clusters before extracting surface- 
brightness and temperature profiles. The aim of this fil- 
tering is to remove high-density, low-temperature features 
which are presumed to be too faint to detect, or perhaps 
even non-physical artifacts of the simulations. We note 
that profiles of unfiltered, individual simulated clusters 
(see, for example, their Figure 2) show significant depar- 
tures from the filtered soft-band surface brightness profile, 
especially at r > 0.5r2oo- The magnitude of these devia- 
tions is at least as great as the north-south difference we 
see in Abell 1795. Similar fluctuations are not seen in the 
hard-band profiles of the simulated clusters. We conclude 
that the north-south surface brightness difference we ob- 
serve is similar to features present in these (unfiltered) 
simulated clusters. 

4-2. Temperature, Density and Entropy Profiles 

The temperature profiles to the north and south of the 
cluster are shown in Figure 10, along with results from 
Chandra and XMM-Newton. We are able to measure the 
cluster temperature past the limit of XMM-Newton and 
Chandra data (about 1.1 Mpc or 15') and find that the 
temperature continues to decline to the largest radius at 
which emission is detectable. At the outermost radial bin, 
spanning 1.3 < r < 1.9 Mpc, the temperature is about one- 
third of the peak cluster temperature observed by Suzaku. 

The radial values and 90% error bars shown in Figure 10 
are estimated from the xissimarf gen ray-tracing simula- 
tions described in Section 2.3.1. These radial bins more 
accurately reflect the real locations of photons originating 
from the cluster, and we use these values for the profile 



fitting described in the remainder of this section. Such 
corrections are less straightforward for the surface bright- 
ness profiles shown in Figures 8 and 9, since the emission 
in that case is from a combination of sources (i.e., cluster 
and cosmic background) that have different spatial distri- 
butions. 

The Suzaku data do not resolve the cool core at r < 2', 
but at larger radii we find good agreement between the 
Suzaku and XMM-Newton in the radial range covered by 
both observatories. Note that here we have used the most 
recent analysis of the XMM-Newton data by Snowdcn 
(private communication), which uses the latest EPIC-pn 
calibrations and provides somewhat lower temperatures 
than those reported earlier (Piffaretti et al. 2005; Snowden 
et al. 2008). The Chandra temperatures reported by 
Vikhlinin et al. (2006), however, are generally higher than 
the Suzaku (and XMM-Newton) temperatures outside the 
cool core, especially at r < 7.5'. The reason for this dis- 
crepancy is not entirely clear. In principle, scattering by 
the telescope of photons from the bright, cool core over the 
central few arcminutes would tend to reduce the apparent 
temperature. One might expect the scattering to be more 
significant for Suzaku, and, to a lesser extent, for XMM- 
Ncwton than for Chandra, given the lattcr's much sharper 
point response function. If this effect were important at 
radii as large asr = 7.5', however, it would be surprising 
if Suzaku and XMM-Newton data agreed as well as they 
do, given their very different point response characteris- 
tics. A more plausible explanation may lie in the known 
discrepancy between Chandra and XMM-Newton temper- 
atures for relatively hot (kT > 4 keV) clusters. 7 Chandra 
reports temperatures that are systematically higher than 
those from XMM-Newton for such clusters, and the mag- 
nitude of the discrepancy at kTxMM ~ 5 keV is very close 
to that shown in Figure 10 at r = 5'. The discrepancy is 
smaller at lower temperatures. Uncertainties in the thick- 
ness of a thin (~ 10-20 A) contamination layer on the 
Chandra mirrors likely account for this discrepancy. We 
note that this effect is unlikely to explain the discrepancy 
at r > 12', where kTxMM ~ 4 keV, a temperature at which 
this calibration error is thought be quite small. We also 
note that the azimuthal range covered by the Chandra 
observations (Vikhlinin et al. 2005) differs from that we 
observed with Suzaku at these large radii. In principle, 
therefore, the discrepancy between Chandra and Suzaku 
results at this large radius could reflect azimuthal tem- 
perature variations in the cluster. In any event, we are 
principally concerned here with the cluster temperature 
profile at r > 5'. Given the good agreement between the 
Suzaku and XMM-Newton temperature profiles in this ra- 
dial range, we base our subsequent analysis on the Suzaku 
temperatures shown in Figure 10. 

Fitting a power law (assuming T cx r~ 7 ) to the pro- 
jected temperature at r > 7.5' yields consistent slopes in 
the north and the south, and 7 = 0.9 ± 0.3 for a joint fit 
to all the data in this radial range. This is similar to the 

7 L. David, 2007, http://cxc.harvard.edu/ccw/proceedings/ 
07_proc/presentations/david/ 
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value 7 = 0.98 ±0.07 reported by George et al. (2008) from 
Suzaku observations of the outskirts of the cluster around 
PKS 0745-191. 

We deprojected the temperature profiles to the north 
and south of the cluster separately using the smaug model 
in XSPEC. The smaug algorithm (Pizzolato et al. 2003) 
fits analytic functions for the deprojected temperature, 
density, and abundance profiles, and compares the repro- 
jected spectra to the data. We fit spectra in the radial 
range 5-17.5' (0.37-1.28 Mpc), although the profiles were 
extended to 2 Mpc for reprojection. Since the emission 
north of 17.5' is evidently not part of a spherically sym- 
metric component, we excluded it from the deprojection. 
We assumed profiles of the form / = / [1 + (r/r c ) 2 ]~ e , 
where /o, r c and e were allowed to vary independently for 
the temperature and density profiles, while fo and e were 
allowed to vary for the abundance profile. The abundance 
profile r c parameter was frozen at 10 kpc, essentially pro- 
ducing a power-law dependence over the fitting region. 
We assumed a mekal model for the plasma emission, and 
included Galactic and extragalactic components identical 
to those described in Section 2. Parameter errors were 
estimated using MCMC techniques identical to those de- 
scribed in Section 2.3.3. 

The best-fit deprojected power-law temperature profiles 
are shown as continuous curves in Figure 11, along with 
the measured (projected) temperatures. The shaded re- 
gions surrounding the curves indicate 90% confidence in- 
tervals. The best-fit model parameters for deprojected 
temperature, hydrogen density, and abundance profiles 
are listed in Table 7. Figure 11 shows that the model de- 
projected temperatures track the observed projected tem- 
peratures quite closely. As expected, the deprojected tem- 
peratures arc slightly higher, but, as suggested by George 
et al. (2008), the difference is relatively small, especially at 
large radius. Note that formally the fit quality is accept- 
able in the south but marginal in the north (see Table 7). 

As expected from the projected temperatures (see 
Figure 10), the deprojected Suzaku temperatures are 
somewhat lower than those inferred from Chandra 
data (Vikhlinin et al. 2006). Moreover, our deprojected 
temperature profile falls more rapidly with increasing ra- 
dius than does the Chandra profile. For example, at 
r = 1.3 Mpc « r5oo, our (average) model profile has a loga- 
rithmic slope -dhxT/dhxr = 7 = 0.77±q Jg while Vikhlinin 
et al. (2006) find 7 = 0.48 ±0.15. (Here we cite 68% con- 
fidence errors.) In interpreting this comparison one must 
again bear in mind that the Suzaku and Chandra observa- 
tions sample different ranges in azimuth around the clus- 
ter. Our measured slope is also steeper than those found 
for simulated clusters by both Hallman et al. (2007), who 
fit a functional form with logarithmic slope 7 ~ 0.4 ±0.2 
at r = rsoo, and by Roncarelli et al. (2006), who report 
7 ~ 0.4-0.6 at r < r2oo- In the latter case the steeper 
slope is associated with simulations which include radia- 
tive cooling. 

Our deprojected density profile approaches a power 
law at large radius, and the results in Table 7 imply 
-d\nn H /d\nr = a = 2.27 ±0.07 at r = 1.3 Mpc. This 
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Fig. 11. Deprojected temperature profiles (solid lines) and 
90% confidence intervals (shaded regions) plotted with the 
observed projected temperatures from Figure 10 (points). 
Results for the north are in black, those for the south are 
in red. Only the annuli between 5—17.5' were included in the 
deprojection. 

is in agreement with the density profile slope reported 
by Vikhlinin et al. (2006) from Chandra data: their fit- 
ting formula is reasonably well-approximated by a power 
law with index a = 2.21 over the radial range we consider. 
Our assumed forms for density and temperature imply 
that the ICM should be polytropic at large radius, with a 
best-fit index of V = 1.3+^. 

Finally, we have computed the implied entropy index, 
(s = kT/rie^ 3 ) as a function of radius, assuming n e = 
1.2 x riff. The result is shown, with 90%-confidence un- 
certainties, in Figure 12. There we also show the entropy 
profile expected from hierarchical structure formation de- 
rived from (non-radiative) simulations by Voit (2005); this 
curve increases as s ~ r . Our average entropy profile in- 
creases more slowly, as s ~ r o.74±o.20^ 

The relatively low temperatures and entropies we ob- 
serve in the outskirts of Abell 1795 may be an indication 
that the plasma we detect is not in hydrostatic equilib- 
rium in the cluster's gravitational potential. We note that 
our Suzaku observations sample two restricted sectors of 
the cluster at this radius, and that data from other ob- 
servatories sample different sectors. If the ICM is not in 
hydrostatic equilibrium at this radius, one might expect 
azimuthal variations in ICM density and temperature. For 
example, as we noted earlier in Section 4.1, the simulations 
of Roncarelli et al. (2006) produce (unfiltered) clusters 
with distinct cool clumps at r < r^oo- These may be the 
result of relatively cool, low-mass substructures falling to- 
ward the cluster center. Moreover, Mahdavi et al. (2008) 
infer modest departures from hydrostatic equilibrium at 
radii r25oo <r < r§QQ from a comparison of cluster masses 
determined from X-ray and wcak-lensing data. This result 
is at least qualitatively in agreement with our interpreta- 
tion of the temperature and entropy profiles of Abell 1795. 
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north 

0.9 
+0.12 



south 



16.2t^ 

n 4o +012 

u - w -o.io 

hO.09 

-0.05 

a - UD -0.09 
1 on+0.03 
1 ' zz -0.03 
+0.01 
-0.01 

-2.91 
+0.06 
-0.12 

0.01 



average 



3Q+ 012 

u - oy -0.09 



kT 



fo (keV) 
c 

r c (Mpc) 



12 

0.38_ c 
19+° 10 

u - la -0.06 



0.13] 



J - uo -0.08 
L17 +0.03 

0.30 



0.16 



+0.10 
0.05 

J - u *-0.08 

1 20+ 003 

x - zu -0.03 
+0.01 
0.01 
+0.24 



/o (10 



-3 



(Mpc) 



0.03 
+0.111 
-0.01 



0.30 



4 Q1+0.2V 
^• Oi -3.02 

S6+ 05 

U - OD -0.14 

0.01 



0.30 



abund fo (solar) 
e 

r c (Mpc) b 



4.40 
0.36 



4.36_ 2 97 
36+ 06 

U - OD -0.13 

0.01 



xVdof 



a Profiles are of the form / = /o[l + (r/r c )' ! ] 
b The abundance core radius was fixed. 
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Fig. 12. Entropy profiles (solid lines) and 90% confidence 
intervals (shaded regions). Results for the north are in black, 
those for the south are in red. Plotted is the entropy index 
s = kTn~ 2/3 in keV cm 2 . Only annuli between 5-17.5' were 
used in the deprojection to determine the profiles. The blue 
line and shaded region are from simulations reported by Voit 
(2005) for non-radiative, hierarchical cluster formation. 

4.3. Cluster Mass 

The X-ray surface brightness and temperature profiles 
can be used to estimate the cluster mass as a function of 
radius if we assume that the ICM is hydrostatic equilib- 
rium. In this case, following Vikhlinin ct al. (2006), 



M{<r) 



-3.68 x 10 13 M ( 



+ 





dlogT 



T(r) 



lkcV 



IMpc 



dlogr dlogr 



(1) 



Here p g is the mass density of the ICM and it is assumed 
that the mean molecular weight of the ICM /1 = 0.62. 
Chandra (Vikhlinin et al. 2006) and XMM-Ncwton (Ikebe 
et al. 2004; Piffaretti et al. 2005) have been used to esti- 
mate mass profiles for Abell 1795 at radii up to r < 1.3 



Mpc. In this section we compare mass constraints ob- 
tained from our Suzaku observations with these measure- 
ments. 

Adopting the mean of the north and south profiles for 
the deprojected temperature and density (see Table 7), 
together with the mean of the north and south tempera- 
tures measured in the final radial bin of each deprojected 
profile, we use Equation 1 to estimate the cluster mass. 
We obtain M = i.ltol * 10 14 M Q within r = 1.27 Mpc, 
the value of r$oo reported by Vikhlinin et al. (2006). The 
68% confidence uncertainties quoted here are dominated 
by uncertainties in the temperature measurement. 

Our mass estimate is somewhat lower than others re- 
ported for Abell 1795. Vikhlinin et al. (2006) find M = 
6.0 ±0.5 x 10 14 M Q at r = 1.27 Mpc, while Piffaretti et al. 

(2005) report M = 5.2 ±0.4 x 10 14 Af Q at r = 1.16 Mpc. 
Formally, these exceed our mass estimate by 2.7 a and 
1.7(7, respectively. The discrepancy is almost entirely due 
to the difference in temperature estimated at these radii. 
For example, the best-fit model used by Vikhlinin et al. 

(2006) for the deprojected temperature predicts T = 3.8 
keV, while we estimate T = 2.9^ keV at r = 1.27 Mpc. 



Our lower mass estimate implies rsoo = 1-08 



+0.07 
^0.03 



Mpc, 



somewhat smaller than the Chandra estimate of rsoo = 
1.27 ±0.04 Mpc (Vikhlinin et al. 2006), and consistent 
with the XMM-Newton value r 500 = 1.16 ±+0.05 Mpc 
(Piffaretti et al. 2005). Similarly, we find r 200 = 1.52±^;^ 
Mpc, smaller than both the Beppo-Sax estimate of r2oo = 
2.14 ±0.46 (Ettori, Dc Grandi & Molendi 2002), and the 
value (f2oo =1-9 Mpc) we estimate from theoretical scal- 
ing relations (Evrard, Metzler & Navarro 1996; Arnaud & 
Evrard 1999). Finally, given that r v ~r W o in our assumed 
cosmology (Bryan & Norman 1998), our mass profile im- 
plies r v = 2.0 Mpc. 

The rapidly falling temperature profile we observe im- 
plies a rapidly falling mass density profile and a slowly 
rising integrated mass profile (sec Equation 1). A power- 
law approximation to our best-fit integrated mass pro- 
file in the region 1 < r < 2 Mpc yields M(< r) = 3.7 x 
10 14 A/ Q ( T j^) - 24 . This is much flatter, for example, 
than the slope expected from an NFW model with scale 
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radius r s = 0.385 Mpc reported by Vikhlinin et al. (2006). 
5. Summary and Conclusions 

We have presented results of Suzaku observations of two 
distinct regions in the outskirts of Abell 1795. We detect 
X-ray emission to r = 1.3 Mpc in both regions, and in 
the north trace the hot ICM to r = 1.9 Mpc > r 20 o- We 
find that the X-ray surface brightness at 1.3 < r < 1.9 
Mpc is significantly higher in the north than in the south. 
We measure the run of temperature with radius at r > 
^2500 an d find that it falls relatively rapidly (Tdeprojected oc. 
r -0 - 9 ) and reaches a value about one-third of its peak at 
the cluster outskirts. 

One possible interpretation of our observations is that 
the ICM is not in hydrostatic equilibrium in the radial 
range r rMQ < r < r 2 oo, and that in the north we are seeing 
infalling plasma at relatively low (kT ~ 2 keV) tempera- 
ture. In support of this picture we note that the X-ray 
surface brightness is neither azimuthally symmetric nor 
(in the north) falling monotonically with increasing ra- 
dius. Moreover, the rapidly falling deprojected temper- 
ature profile would require that there is relatively little 
gravitating mass at r > 1.3 Mpc if the ICM is polytropic 
and in hydrostatic equilibrium, a result that conflicts 
with extrapolations of the best-fit NFW profiles derived 
from high-resolution Chandra observations (Vikhlinin et 
al. 2006). 

Simulations (Roncarclli et al. 2006) predict, and re- 
cent lensing and X-ray observations (Mahdavi et al. 
2008; George et al. 2008) provide evidence for, modest 
deviations from hydrostatic equilibrium in the ICM at 
r 2500 < r < ?*2oo- Our observations of Abell 1795 suggest 
that in some cases these deviations may be quite signif- 
icant even in a cluster which appears to be relaxed at 
smaller radii. Additional observations of a variety of clus- 
ters, with more complete azimuthal coverage than we have 
obtained, will improve understanding of the state of the 
ICM in cluster outskirts, and therefore of cluster masses, 
structure and formation mechanisms. 

Finally, we note that these and other observations of 
the low-surface-brightness outskirts of clusters have been 
made possible by Suzaku's relatively low and stable non- 
X-ray background. The surface brightness sensitivity of 
our observations, -Bo.5-2fceV > 1-8 x 10 -12 erg s _1 cm~ 2 
deg~ 2 (3 a), is less than 20% of the total cosmic back- 
ground in this spectral band. Our sensitivity was ul- 
timately limited by time-variable solar wind charge ex- 
change emission and by Poisson variations in the extra- 
galactic source density. We demonstrated here that vari- 
ations in the geocoronal contribution can be modeled to 
this level from simultaneous observations of cluster-free 
background regions. To reduce background source den- 
sity variations below this level in the relatively small clus- 
ter regions (0.01-0.03 deg 2 ) we required information from 
the XMM-Newton source catalog. Future Suzaku obser- 
vations of cluster outskirts would benefit from coverage of 
larger (~0.1 deg 2 ) solid angle in each radial bin, especially 
if deep X-ray source catalogs are not available. 
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